Vibrational origin of the fast relaxation processes in molecular glass-formers 
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We study the interaction of the relaxation processes with the density fluctuations by molecular 
dynamics simulation of a flexible molecule model for o-terphenyl (oTP) in the liquid and supercooled 
phases. We find evidence, besides the structural relaxation, of a secondary vibrational relaxation 
whose characteristic time, few ps, is slightly temperature dependent. This i) confirms the result by 
Monaco et al. [Phys. Rev, E62, 7595 (2000)] of the vibrational nature of the fast relaxation observed 
in Brillouin Light Scattering (BLS) experiments in oTP; and ii) poses a caveat on the interpretation 
of the BLS spectra of molecular systems in terms of a purely center of mass dynamics. 

PACS number(s): 64.70.Pf, 71.15.Pd, 61.25.Em, 61.20.-p 
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After the development of the Mode Coupling Theory 
(MCT) §, aimed to give a microscopic interpretation to 
the slowing down of the dynamics which takes place at 
the liquid-to-glass transformation, different experimen- 
tal and numerical studies have been devoted to test its 
predictions, especially through the study of the density 
fluctuations correlation function, or of its space Fourier 
transform F(Q, t). These studies can be assigned to two 
broad classes, according to the value of the investigated 
wavevector Q. On one side the numerical techniques and 
the inelastic neutron scattering spectroscopies, that ac- 
cess the high Q region, give a picture where the MCT 
predictions are basically demonstrated. On the other 
side the low Q techniques as, for example, Brillouin Light 
Scattering (BLS), depolarized light scattering, or dielec- 
tric spectroscopy, depict a much less clear situation. For 
example, the recent observation of a flat background 
(constant loss) in the dielectric absorption spectra at few 
GHz ||], seems to contradict the MCT predictions on 
the same systems where the high Q techniques confirm 
these predictions. A similar conflictual situation is also 
found in the case of the BLS experiments ||, in par- 
ticular as far as the so called /3-region of the MCT is 
concerned. This is the time region where the first de- 
cay of F(Q,t), that could be ascribed to the dephasing 
of the microscopic vibrational dynamics, merge with the 
earliest part of the decay associated to the structural (a) 
relaxation process. The MCT makes specific predictions 
on the shape of the F(Q, t) in this region, but these pre- 
dictions are not in agreement with the BLS, low Q, mea- 
surements 0] whereas the high Q neutron spectroscopy 
H and molecular dynamics data support the validity 
of the MCT. 

Recently, studying the BLS spectra of o-terphenyl 
(oTP), one of the prototypical fragile glass formers, some 
of us found evidences -besides the usual structural or a 
relaxation process- of a secondary relaxation, with a char- 
acteristic time Tf lying in the ten picosecond time-scale 



and weakly temperature dependent. This process lies in 
the frequency region where the features associated to the 
(3 region are expected to show up. By comparing the 
BLS spectra in the glassy and in the ordered crystalline 
phases, and by comparing the effect of this relaxation 
process on the longitudinal and transverse sound waves, 
it has been suggested the vibrational nature of the ob- 
served relaxation ||] . If confirmed, this suggestion would 
indicate the reason why the MCT predictions are not 
always retrieved in the light scattering and dielectric ex- 
periments. Indeed the MCT is, up to now, formulated 
for monatomic or rigid molecules, the complexity of real 
non-rigid molecules -with all the processes associated to 
the internal vibrational dynamics- not being taken into 
account. The presence of a vibrational relaxation would 
indicate that the molecular glass formers are not suitable 
benchmarks to test the MCT on the short time scale. In- 
deed the presence of a vibrational relaxation dynamics 
lying in the 1-100 ps time-scale would mask the spectral 
features predicted in the /3-region. This considerations 
could explain the lack of consistency between the MCT 
predictions and the low Q experiments. 

In this letter we present a molecular dynamics simu- 
lation investigation of the relaxation processes active in 
oTP in a wide temperature range, spanning the normal 
and supercooled liquid region. Using the recent flexi- 
ble molecular model for oTP ||[l0]], we investigate the 
coupling between the density fluctuations and the intra- 
molecular (vibrational) degrees of freedom. We find that 
i) there is a strong coupling between the density fluctu- 
ations and the vibrational excitations, and ii) that this 
coupling gives rise to a vibrational relaxation in the 10~ n 
s time-scale. 

These findings strongly support the hypothesis sug- 
gested in Ref. ^ and confirm the existence of a vibra- 
tional relaxation in oTP in the same frequency region 
where the signatures of the MCT /3-region are expected. 
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FIG. 1. Time dependence of the total intra- molecular po- 
tential energy after the density jump starting from a well 
equilibrated configuration at T = 325 K. The data have been 
obtained averaging over 2,000 realizations and the error bars 
represent the ±cr standard deviation of the different runs data. 

Moreover, by analyzing the effects of an external per- 
turbation on the individual contributions to the intra- 
molecular interaction potential, we identify the phenyl- 
phenyl stretching as the internal vibration mainly respon- 
sible for such a relaxation. 

The absolute value of the relaxation time Tf turns out 
to be a factor w 3 shorter than the experimental value, 
a discrepancy that is tentatively explained on the basis 
of the classical nature of the MD simulations. Finally, 
we investigate the temperature dependence of Tf, find- 
ing a behavior very similar to that observed in the BLS 
experiment 1 1 1 . 

We have studied a microcanonical system composed 
by 32 molecules (96 rings, 576 Lennard Jones interaction 
sites) enclosed in a cubic box with periodic boundary 
conditions. In the utilized model [p| the oTP molecule 
is constituted by three rigid hexagons (phenyl rings) of 
side L a = 0.139 run. Two adjacent vertices of the parent 
ring are bonded to one vertex of the two lateral rings by 
bonds of equilibrium length = 0.15 nm. Each vertex 
of the hexagons is occupied by a fictious atom of mass 
Mch — 13 a.m.u. representing a carbon-hydrogen pair 
(C-H). In the isolated molecule equilibrium position the 
lateral rings lie in planes that form an angle of about 54° 
with respect to the parent ring's plane. The interaction 
among the rings pertaining to all the molecules is de- 
scribed by a site-site pairwise additive Lennard- Jones 6- 
12 potential, each site corresponding to the six hexagons 
vertices. Moreover, the three rings of a given molecule 
interact among themselves by an intra-molecular poten- 
tial, such potential being chosen in such a way to pre- 
serve the molecule from "dissociation", to give the cor- 
rect relative equilibrium positions for the three phenyl 
rings, and to represent at best the isolated molecule vi- 
brational spectrum. In particular it is written in the form 
^intra = J2k c k^k where each term 14 controls a partic- 
ular degree of freedom. The contributions to the intra- 
molecular potential relevant to the present study are: i) 



the stretching along the central ring - side ring bonds (S), 
ii) the bending of the central ring - side rings bonds {B), 
and iii) the in-phase {R\) and out-of-phase (R2) rotation 
of the lateral rings along the ring-ring bond p2| . The 
first two terms are modeled by springs, and the "inter- 
action" with the other degrees of freedom (anharmonic- 
ity) is introduced by the site-site LJ potential of differ- 
ent molecules. The third term, on the contrary, has been 
modeled in more realistic way. In this case the relevant 
variables are the two angles { ( &i, ( l > 2} between the nor- 
mals to the lateral rings and the parent ring plane. An 
ab-initio calculation of the single molecule potential en- 
ergy surface as a function of these two angles has shown 
that exist two iso-energetic configurations separated by 
an energy barrier of height V s /kB — 580 K. The nature 
of the rotational motion at the temperatures of interest 
can be summarized as follows: the two side rings can 
pivot in phase around the bonds crossing from one mini- 
mum to the other degenerate one; at the same time they 
can perform librational out-of-phase motions. In order to 
represent this potential surface we express the in-phase 
rotation of the two side rings with a high-order (6th) 
polynomial and the out-of-phase rotation by a quadratic 
(harmonic) potential energy. Other details of the intra- 
molecular and inter-molecular interaction potentials, to- 
gether with the values of the involved constants, are re- 
ported in Ref. Previous studies on the temperature 
dependence of the self diffusion coefficient jj| and of the 
structural (a) relaxation times PJlO| indicate that this 
molecular model is capable to quantitative reproduce the 
dynamical behavior of the real system, but the actual 
simulated temperatures must be shifted by ~ 20 K up- 
ward. In what follows, as our aim is to compare the sim- 
ulation results with the experiments, the reported MD 
temperatures are always shifted by such an amount. 

To investigate the effect of a long wavelength density 
fluctuation, as those probed in a BLS experiment, on 
the intra-molecular vibrational dynamics we proceed as 
follows. After an equilibration run at a given tempera- 
ture, we make a sudden density variation of the system 
("crunch"), then we follow the time evolution of the in- 
tramolecular potential energy during the subsequent evo- 
lution. As an example, in Fig. ^ we report such a time 
evolution for T=325 K and averaged over 2,000 statisti- 
cally independent crunches. The value of the internal po- 
tential energy is close to that pertaining to 12 harmonic 
oscillators (the number of internal degrees of freedom), 
6i?T=16.2 kJ/mol, the slight deviation being associated 
to the anharmonicity present in the parameterizations of 
some of the internal degrees of freedoms. It is clear that, 
after the density change, the vibrational energy relaxes 
toward its new equilibrium value in a time-scale of «6 
ps, a value not far from the experimental one: Tf «20 
ps. We do not expect a much better agreement between 
these two values as i) the intra-molecular potential model 
has been parameterized to represent the vibrational 
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FIG. 2. (A) Signals (main panel) and relaxation function 
(inset) of the potential controlling the on-phase pivoting of 
the lateral rings at T = 370 K and A = 1.2. (B) As above for 
the out-phase librations. (C) As above for the bending of the 
central ring - side rings bond. (D) As above for the stretching 
between the sides and central ring along the molecular bonds; 
it is worth noting that V s decays on a time scale much longer 
than that of the other terms shown above. 

spectrum of the isolated molecules -and not on the dy- 
namics in the condensed phase-, and ii) the intramolecu- 
lar vibrational frequencies have values up to «500 cm -1 
(«720 K) and therefore, at the investigated temperature, 
the quantum effects (population of the vibrational levels) 
certainly play a relevant role, and they are not considered 
in the present classical simulation. 

Overall, we conclude from Fig.[l] that there is defini- 
tively a coupling between the density fluctuations and 
the vibrational degrees of freedom -i. e. a vibrational re- 
laxation is active in the system- and that the relaxation 
time favorably compares with the experimental findings. 

In order to identify which one of the internal degree 
of freedom(s) is more efficiently coupled to the density 
fluctuations we cannot use the crunch technique, as a 
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FIG. 3. Relaxation functions Rs(t) of the bond stretching 
potential at the four indicated temperatures. 



density jump simultaneously affects all the intra- 
molecular degrees of freedom. In order to selectively 
perturb a specific vibration, we proceed as follows. Af- 
ter the equilibration run at a given temperature, the 
value of one of the elastic constants Cfc is sudden scaled 
(c fe — > c* k = Ac fc ) [l3|Jl^1 with A positive and small |]|. 
After this perturbation is applied, we follow the time evo- 
lution of the specific term of the intra-molecular interac- 
tion potential energy and we measure the time needed to 
the energy to relax towards its new equilibrium value. As 
an example, in Fig ^ (A-D), we report the time evolution 
of some of the terms Vt (k = Ri, R2, B, S) for the se- 
lected temperature T=370 K (open squares). The value 
of the energies measured in runs without the perturbation 
(open triangles) are also shown for comparison. Finally, 
in the inset, the difference Rk(t) (normalized to unity at 
t=0) between the perturbed and unperturbed time evo- 
lution are shown. As can be seen, all the potential term 
but Vs decay on a very fast time-scale (wlOO fs or less fol- 
lowed by an oscillating tail that vanish in «0.5 ps). The 
term Vs on the contrary decays on a much longer (> 1 
ps) time scale. This term, which controls the stretching 
along the two bonds connecting the parent ring with the 
two side rings, is therefore the main candidate for the 
vibrational relaxation observed at few ten ps. It is worth 
noting that, at variance with the other cases, the per- 
turbed signal does not relax at the same value pertaining 
to the equilibrium signal; this indicates the presence of 
a relevant anharmonicity affecting this degree of freedom 
and that its equilibrium value depends on the density. 

The functions Rs(t) are reported in Fig. || at selected 
temperatures, together with their best fit to an exponen- 
tial decay. A slight but evident temperature dependence 
of the relaxation time is present. The T dependence of 
the relaxation time derived from these fits is reported in 
Fig [|, together with the relaxation times for the fast 
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FIG. 4. Experimental (open circles, From Ref. [11]) and 
MD (full diamonds) results for r/(T) reported in an Arrhc- 
nius plot. The MD data have been multiplied by a common 
factor F = 6 in order to show the T behaviour similar to that 
of the experimental data. 



process experimentally determined in Ref. 11 1. In this 
figure the MD data have been multiplied by a common 
factor F = 6 in order to emphasize that their T behaviour 
is very similar to that of the experimental data. Both 
sets of data show an Arrhenius behaviour at low T (with 
activation energy Ai?=0.28±0.01 kJ/mole) and a differ- 
ent (steepest) trend at high T, where the fast relaxation 
process merges with the structural one. The agreement 
between the T behaviour of the two sets of data is re- 
markable. On the contrary, there is a discrepancy of a 
factor of about 6 between the MD and the experimental 
time-scales. This difference could be explained, as no- 
ticed above, considering a non perfect parameterizations 
of the intra-molecular interaction potential and the fact 
that we expect strong quantum effects on the studied 
process. Moreover, in the present case, further differ- 
ences can arise from the different " perturbation" utilized 
in the MD simulation (changes of force constants) and in 
the experiment (density fluctuations). 

In conclusion, by studying the coupling of external per- 
turbations (a density fluctuation and a fictious coupling 
constant modification) with the intra-molecular vibra- 
tions in a flexible model of oTP molecule we demonstrate 
i) the existence of a vibrational relaxation process with a 
non negligible strength and a relaxation time in the few 
ps time-scale; ii) that this relaxation process is mainly as- 
sociated to the phenyl-phenyl stretching; and iii) that the 
relaxation time shows a temperature behavior very sim- 
ilar to that of the experimentally determined Tf [MHO . 
The quantitative discrepancies between the simulate and 
the experimental relaxation times can be tentatively as- 
signed to the quantum nature of the real vibration, for 
which, at room temperature, hu> v /kBT w 2. The present 



findings give strong support to the vibrational origin of 
the fast relaxation process observed in oTP. It is tempt- 
ing to generalize this conclusion to other systems where 
such a fast process has been observed (PC jl6| and PB 
p^t). Giving the depicted scenario, one should be ex- 
tremely careful in drawing conclusion on either the va- 
lidity, or the failure, of the MCT by the analysis of the 
light scattering or dielectric spectra. This is especially 
true when analyzing the MCT /3-region, as it lies at fre- 
quencies that coincide with the typical fast process ones 
and as its main features (susceptibility minima, a and b 
exponents, knees, etc.) can be masked by the fast process 
itself. 
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